Abstract. We report the results of a study of the ra- 
dio surroundings of the superluminal microquasar GRO 
J1655— 40. We have searched for extended continuum 
structures that might be indicative of the presence of a 
supernova remnant (SNR) associated with the formation 
of the compact object in the binary system. We also car- 
ried out Hl-line observations of the region looking for a 
local minimum created by an explosive event. Our results 
indicate that there is, in fact, a bubble in the large-scale 
" l ' HI distribution around GRO J1655— 40. We suggest that 

. this structure might be created by the original supernova 
<^ ' explosion occurred a few hundred thousand years ago and 

, whose signatures can be traced by the overabundance of 
5h ' a-elements recently found by Israelian et al. (1999) in the 

, companion star. 
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1. Introduction 

The soft X-ray transient GRO J1655— 40 is one of the 
best black hole candidates in the Galaxy (Mirabel & Ro- 
driguez 1999). It is a low mass X-ray binary whose com- 
pact object has a mass in the range 5.5 — 7.9 M Q (Shah- 
baz et al. 1999). The system is located at 3.2 kpc from 
the Earth and presents frequent superluminal ejections at 
radio wavelengths (e.g. Tingay et al. 1995, Hjellming & 
Rupen 1995). 

Recently, Israelian et al. (1999) found convincing ev- 
idence of a supernova origin for the black hole in GRO 
J1655— 40. They detected a clear overabundance of a- 
elements in the secondary star, whose internal tempera- 
tures are not high enough as to synthesize them. These el- 
ements seem to have been produced in a massive compan- 
ion (with 25 — 40 Mq) which exploded both contaminating 
the secondary star with matter enriched with nucleosyn- 
thetic products and creating the black hole. Israelian et al. 
(1999) estimated that the explosion should have occurred 
less than 10 6 years ago. 

In this Letter, we present the results of a study of the 
radio environment of GRO J1655— 40. Our aim is to detect 
any existent signature of an explosive event in this region. 
A supernova explosion should have introduced modifica- 
tions in the ISM around the binary system, displacing ma- 
terial and creating a low density bubble. Additionally, the 
expanding shock front should accelerate electrons up to 
relativistic energies, producing a shell-type source in the 
radio continuum. If the age of the system is less than 10 6 
years, these observational features might be detectable. 

2. Data analysis and new observations 

2.1. Radio continuum 

We have used data from the 4.85-GHz PMN survey (Con- 
don et al. 1993) in order to produce a large-scale map 
(~ 1° x 1°) of the surroundings of GRO J1655-40. Al- 
though these data are optimized for sizescales of less than 
40 arcminutes, the good resolution and sensitivity of the 
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Fig. 1. Upper panel: Filtered radio emission at 4.85 GHz 
of the region surrounding GRO J1655— 40. The position 
of the microquasar is marked by a cross. Contours are 
labeled in steps of 10 mJy beam -1 , starting from 35 mJy 
beam -1 . Lower panel: Small-scale, higher resolution, VLA 
image at 1.4 GHz for the region close to GRO J1655— 40. 
Radio contours are in steps of 2 mJy beam -1 starting from 
1 mJy beam -1 . The arrows indicate the directions of the 
observed superluminal ejections. 

images render them a valuable tool to search for SNR 
candidates, especially if they are not of large angular size 
(Duncan et al. 1997). For this work, we apply to the 4.85- 
GHz data an additional filtering process in order to remove 
the galactic diffuse emission on scales larger than 30 ar- 
cminutes (see Combi et al. 1998 for details of the Gaussian 
filtering method). The resulting map is shown in Figure 
1, along with a 1.4-GHz image of the small-scale emission 
of the inner region obtained from the NVSS Sky Survey 
(Condon et al. 1998) at a better resolution. These lat- 
ter observations were obtained with the Very Large Array 
(VLA) in the compact D and DnC configurations. 

2.2. HI observations 

We have performed HI observations towards the area of 
interest with a 30-m single dish telescope located at the In- 
stituto Argentino de Radioastronomia (IAR), Villa Elisa, 
Argentina. The observations were carried out during four 
consecutive sessions on November 20-23, 2000. The re- 
ceiver is a helium-cooled HEMT amplifier with a 1008- 
channel autocorrelator at the backend. System parameters 
and additional details of the observational technique can 
be found in Combi et al. (1998). The HI line was observed 
in hybrid total power mode and the sky was sampled on 
a 0.35° rectangular grid. Each grid position was observed 
during 60 s with a velocity resolution of ~ 1 km s -1 and a 
coverage of ±450 km s -1 . A set of HI brightness temper- 
ature maps (AT rms ~ 0.2 K) were made for the velocity 
interval Av =(—80 km s -1 , +20 km s -1 ). Those maps for 
the interval —42 to —30 km s -1 , corresponding to the HI 
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Fig. 2. HI brightness temperature channel maps (contour 
labels in K) obtained for the velocity range from —42 to 
—30 km s _1 around the position of GRO J1655— 40 (in- 
dicated with a cross). Due to the differential galactic ro- 
tation the expected velocity of the interstellar medium at 
the distance of 3.2 kpc of the X-ray binary is ~ 38 km s _1 



distribution at a distance centered at ~ 3.5 kpc, are shown 
in Figure 2. 

3. Main results 

Figure 1 (upper panel) shows the radio continuum image 
at 4.85 GHz of the surroundings of GRO J1655-40. A 
small (~ 24' x 24') semi-circular radio source, which very 
much resembles a SNR, can be clearly seen in the center of 
the image. This extended, partial shell source is centered 
at (I, b) « (345.14°, +2.48°). It has an integrated flux den- 
sity of ~ 0.5 Jy at 4.85 GHz, although uncertainties due 
to the background filtering make difficult to estimate a 
realistic value. The microquasar is located at a distance 
of ~ 0.15 deg from the center of the radio source. 

Figure 1 (lower panel) shows the VLA image at 
1.4 GHz (Condon et al. 1998) of the vicinity of GRO 
J1655— 40. The microquasar was undetectable at the 
epoch of the observations and its position is marked with a 
cross. There are seven point sources near the microquasar 
with fluxes above 5 mJy. They have been labeled from 
SI to S7 in the image (notice that the SNR candidate is 
not visible in this map because these data are optimized 
for point source detection). The bright, compact source 
located at (I, b) ~ 345.1°, +2.37° is identified in the PMN 
catalogue as PMN J 1654-3949. The sources SI to S6 have 
no entry in any point source catalog at present. The mea- 
sured characteristics of all these sources are summarized 
in Table 1. 

Figure 2 clearly shows the development of a cavity 
around GRO J1655-40 in the HI maps from -42 to -30 
km s _1 . Although the cavity do not form a shell due to 
the steep density gradient towards the galactic plane, the 
existence of a local minimum in the HI distribution at 
the position of the microquasar is evident. This kind of 
"horse-shoe" morphology has been observed around early- 
type stars with very strong stellar winds (e.g. Benaglia & 
Cappa 1999), where a strong energy release in a highly 
inhomogeneous medium has occurred. 

In Figure 3 we present the integrated column density 
map for the velocity interval from —40 to —30 km s _1 . 
Standard galactic rotation models (Fich et al. 1989) indi- 
cate that this HI hole is located at a kinematic distance of 
3 to 4 kpc, in good accordance with the estimated distance 
to GRO J1655-40. 



Table 1. Point radio sources near GRO J1655— 40 



Source 


(l,b) 


F 1A GHz ID 




(deg, deg) 


(mJy) 


SI 


(344.87°, +2.62°) 


24.9 


S2 


(344.92°, +2.54°) 


12.8 


S3 


(344.96°, +2.46°) 


5.43 


S4 


(344.98°, +2.42°) 


4.48 


S5 


(345.02°, +2.38°) 


37.7 


S6 


(345.03°, +2.49°) 


5.83 


S7 


(345.10°, +2.37°) 


10.6 PMN J1654-3949 



Fig. 3. Integrated column density map for the velocity 
range from —40 to —30 km s _1 . Contour labels are in 
units of 10 19 cm" 2 . Microquasar position indicated with 
a cross. 



4. Discussion 

If we assume that the extended radio structure shown in 
Fig. 1, upper panel, is a supernova remnant located at 
the same distance than the microquasar, we get that its 
radius should be ~ 15 pc. The surface brightness of the 
source results to be abnormally low for such a remnant 
(£ 4 .8 GHz ~ 1.3 x 10~ 22 W m- 2 Hz" 1 sr- 1 ). On the other 
hand, the microquasar is displaced about 0.15° from the 
center of the semi-shell source. Although the transverse 
velocity of GRO J1655— 40 is not known, it is constrained 
to be less than 300 km s" 1 (Brandt et al. 1995). At this 
velocity, it would have taken more than 2.7 x 10 4 yr to 
reach its current position. And for a velocity similar to 
the observed radial velocity (~ 114 km s _1 ), the age of 
the system would be <~ 82 500 yr. These timescales are 
incompatible with the size of the remnant and our knowl- 
edge of the ambient density in this region (see below), and 
lead us, consequently, to discard the extended radio source 
in Fig. 1 as possibly associated with GRO J1655— 40. This 
source could be a background SNR or, perhaps, an artifact 
formed by several close and unresolved point sources. 

At larger sizescales there are not other clear contin- 
uum structures that could be associated with a SNR. We 
have inspected larger maps, of ~ 6° x 6°, in search of 
additional signatures, but nothing was found. This is a 
very confused region located nearly towards the Galac- 
tic Center and the detection of weak continuum features 
with low surface brightness is extremely difficult. Even at 
small scales there is a significant number of sources, as it is 
shown in Fig. 1, lower panel. Here, the two closest sources 
to the microquasar, S3 and S6, seem to be at approximate 
equal angular distances, in opposite directions, of the X- 
ray binary. However, it is unlikely that they could be old 
emission knots ejected by the central object, because they 
are on an axis displaced 30° from the axis determined by 
the known ejections (Hellming & Rupen 1995), which is 
also shown in the figure. 
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Contrary to the often misleading continuum data, the 
HI observations can provide more clear signatures of an 
explosive event in this region, because the differential ro- 
tation motion of the material allows a glimpse of the mat- 
ter distribution along the line of sight. The channel maps 
in Fig. 2 show the development in the velocity distribu- 
tion of a cavity in the HI with a "horse-shoe" shape and 
an average radius of ~ 75 pc at a distance of 3.2 kpc. The 
integrated column densities shown in Fig. 3 can be used to 
estimate the mass displaced to form the cavity. Assuming 
that the background density is given by the density at the 
minimum, we estimate that at least 25 000 M Q have been 
removed from this region. Since the average expansion ve- 
locity of the shell is ~ 15 km s _1 , we obtain that a lower 
limit to the original energy release which started the mass 
motion at the center of the cavity is ~ 6 x 10 49 erg. The 
expansion velocity should have been considerably higher 
in the past, in such a way that the real value is probably in 
the range 10 50 — 10 51 erg, what is quite consistent with the 
supernova hypothesis. This is reinforced by the fact that 
there are not known early-type stars with stellar winds 
strong enough to create such a bubble in this sector of the 
sky. It is noteworthy that even higher values for the initial 
energy release are possible if the explosion was caused by 
a hypernovae event as recently proposed by Brown et al. 
2000. 
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The mean number density in this region can be ob- 
tained averaging the column densities within the bubble. 
This provides a rough estimate of the average density of 
the ISM before the explosion. We obtain a value n ~ 0.4 
cm~ 3 . This density and the size of the cavity (R ~ 75 pc) 
indicate that the SNR should be well within the radiative 
phase. An exact estimate of the time elapsed since the ex- 
plosion is not possible without a knowledge of the detailed 
structure of the ISM in this region, but simulations with 
an original energy release of few times 10 50 erg and some 
reasonable conditions for the cloud properties suggest ages 
in the range 10 5 - 10 5 5 yr (Cowie et al. 1981). The con- 
tinuum emission of the remnant, with sizescales of several 
degrees, should have been removed during the filtering of 
the background radiation in the PMN survey. New ob- 
servations with good sensitivity and single dish telescopes 
over a large area are necessary to detect the SNR radio 
structure and help to make a more accurate estimate of 
its age, if the remnant is not already dissipated. 



In summary, we have presented radio continuum and 
HI maps of the surroundings of the black hole candidate 
GRO J1655— 40. A bubble centered at the microquasar 
has been detected in our HI data. The energy necessary to 
create this bubble and its approximate age are consistent 
with the hypothesis that it is the trace left in the ISM by 
the supernova explosion occurred when the black hole was 
formed. 
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